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ABSTRACT 

Observations at submillimeter wavelengths have revealed a population of 
sources thought to be at relatively large redshifts. The position of the 850 
fxm passband on the Ray leigh- Jeans portion of the Planck function leads to a 
maximum redshift estimate of z ~4.5 since sources will not retain their redshift 
independent brightness close to the peak of the Planck function and thus drop out 
of surveys. Here we review evidence that ice absorption is present in the spectra 
of local ultraluminous infrared galaxies which are often taken as analogs for the 
850 fim source population. We consider the implication of this absorption for ice 
induced spectral structure at far infrared wavelengths and present marginal as- 
tronomical evidence that amorphous ice may have a feature similar to crystalline 
ice near 150 /im. Recent corroborative laboratory evidence is supportive of this 
conclusion. It is argued that early metal enrichment by pair instability SN may 
lead to a high ice content relative to refractory dust at high redshift and a fairly 
robust detection of ice emission in a z = 6.42 quasar is presented. It is further 
shown that ice emission is needed to understand the 450 fim sources observed in 
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the GOODS-N field. We are thus encouraged to apply far infrared ice emission 
models to the available observations of HDF 850.1, the brightest submillimeter 
source in the Hubble Deep Field. We suggest that a redshift as large as 13 may 
need to be considered for this source, nearly a factor of three above the usual top 
estimate. Inclusion of the possibility of far infrared ice emission in the spectral 
energy distributions of model sources generally broadens the range of redshifts to 
be considered for submillimeter sources compared to models without ice emission. 

Subject headings: galaxies: high redshift, formation — infrared: galaxies — 
galaxies: individual: IRAS 00188-0856, IRAS 14348-1447, QSO J1148+5251, 
HDF 850.1 



Introduction 



Sources discovered in surveys undertaken using 450/850 fim SCUBA camera mounted 
on the JCMT have a high numb er density per unit solid angle (~ 10"^ deg~^; Sgso > 1 niJy; 
Sanders! I2OOOI : IWang et al.ll2004l ) and are thought to trace the dust emission from sources at 
a range of redshifts. Sources brighter than 1 mJy which have Arp 220-like far infrared (FIR) 
spectral energy distributions (SEDs) are ultraluminous in the /iioo = 0.71, Qm = 0.27, and 

= 0.73 cosmol ogy consistent with cosmic microwave background (CM B) and Type la 
SNe observations (ISpergel et al.l l2003l : iNolta et al.l |2004J : iRiess et al.ll2004l ) unless they are 
at redshifts < 1 so that, thus far, submillimeter surveys do not detect "normal" galaxies at 
high redshift. This should be the case even for gravitational lens amplification up to a factor 
of 4 for 2 mJy limit surveys. 

It is therefore a sensib le approach to consider the properties of local ultraluminous 
infrared galaxies (ULIRGs; ISanders fc Mirabell Il996l ) as a guide to understanding sources 
detected in submillimeter surveys with the proviso that the high redshift universe may not 
conserve these properties. It has already been noted that the number density of ULIRGs 
increases with re dshift at a rate much larger than can be explained using geometric con- 
siderations alone (IKim fc Sanders! Il998l ) and deep field observations using the Hubble Space 
Telescope have found quite obvious morpholo gical evolution at higher redshifts compared to 



local conditions {e.g. !Brinchmann et al 



1998! ). While these particular evidences of evolution 



as a function of redshift do not directly speak to ULIRG SED evolution, they do suggest 
that this proviso cannot be lightly discounted. 



Ice is present in the dense ISM of the Galaxy. It is observed to occur above a threshold 
depth of Ay ~ 3 within molecular clouds but is absent on the cloud edges and in the diffuse 
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(atomic/ionized) ISM (IWhittet et al.lll988[ ). Because ice is not exposed to star light, far 
infrared (FIR) emission from ice is not expected from quiescent clouds. The FIR emission 
from these ice bearing clouds arises from the first few ice- free Ay. To observe ice emission 
in the Galaxy, centrally heated sources are needed. If mid-infrared (MIR) radiation can 
illuminate the ice bearing grains so that they may be heated, then the FIR emission will be 
dominated by ice if the ice mantles provide a moderate fraction of the total grain volume. 
This FIR dominance is owing to gre ater FIR opacity o f ice relative to silicate or carbonaceous 
grains. In the models presented by lAannestadl (119751 ). a mantle-to-core volume ratio of 0.33 
doubles the 60 /um opacity: a 10% increase in grain radius doubles the opacity! 

Protostars in the Galaxy w hich are centrally h eated turn out to be good sources for 



observing FIR ice emission {e.g. iDartois et al.lll998l ). And, there is mounting evidence that 



a significant fraction of local ULIRGs are also centrally heated (Dudlev fc Wvnn- William^ 



1997 



2006a 



Soifer et allbood: llmanishi fc Dudlevlboool : Emanishi fc Malonevll2003l : llmanishi et al.l 



Spoon et al. 20061 ) so that the FIR emission is a result of MIR heating. Thus, if ice 



is present in the FIR emitting regions, its effects on the FIR opacity may be important. 
This situation differs from what occurs in lower luminosity starburst galaxies where the FIR 
emission arises from the surfaces of molecular clouds which are heated by optical and far 
UV rather than MIR light. The interiors of the clouds, where ice presumably exists, are not 
strongly heated by the mainly FIR emission arising from their surfaces. 

If submillimeter galaxies are similar to ULIRGs in more than just their luminosity, 
then they may also sometimes have FIR emission that is powered by MIR light. For redshift 
estimation the question of the similarity between submillimeter galaxies and local ULIRGs is 
as important as the differences that distinguish ULIRGs from their lower luminosity starburst 
counterparts. 

Attempts to estimate redshifts for the submillimeter sources based on their radio fluxes 
coul d be guided by local ULIRG samples. For the 13 ULIRGs brighter than 5.24 Jy at 60 
//m (jSoifer et al.l 1 198 71 : ISanders et al.lll995l ) available to the VLA, the m ean logarithm of the 
ratio of FIR-to-20 cm radio flux is 0.22 dex larger than g=2.43±0.19 JCondon et al.lll99lh 



for sta rbur st galaxies when calculated using radio measurements reported by ICondon et al. 



(1991 


) and 


Yun et al. 


(2001) 



g is a factor 1.6 (in dex) lar ger than for t he low er luminosity star forming galaxies. Using a 



different luminosity binning, iHelou et al.l (119851 ) report a similar trend in increased dispersion 



with increasing luminosity. The implication for redshift estimation is that the range of likely 
redshifts would increase when using ULIRG properties rather than starburst properties. 
And, a trend away from what is seen in local ULIRGs, say through more powerful starbursts 
at earlier times, could have complex implications for estimating redshifts of submillimeter 
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galaxies using radio observations. 

However, broad limits on the redshift distribution derive from submillimeter observa- 
tions alone where the relative infrequency of detecting 450 ^m emission leads to (weak) lower 
limits on the redshifts while the detections themselves at 850 /um lead to upper limits since, 
baring extraordinary luminosities, the turnover from the Rayleigh- Jeans-like portion of the 
SED at high redshift leads to source fluxes dropping below survey detection limits. 

Given the large potential effect of ice on the FIR opacity of submillimeter galaxies, we 
first review the evidence for the presence of ice in centrally heated ULIRGs and estimate an 
expected FIR SED for a case where FIR ice opacity can be expected to be dominant based 
on a rough abundance analysis. We then search for FIR ice spectral features in the available 
data on ULIRGs. We discuss why a high relative abundance of ice during the epoch of 
reionization might be expected and present somewhat firmer evidence for FIR ice emission 
in the z = 6.42 source QSO J1148-I-5251 than can yet be demonstrated for ULIRGs. We 
then take the submillimeter galaxy HDF 850.1 as an example to show how FIR ice emission 
allows a very large redshift estimate. 



2. Ice in ULIRGs 
2.1. Absorption 

Observations of strong absorption features from refractory dust in ULIRGs trace dust 
that is cooler than the temperature of a blackbody emitting at the wavelength of the feature: 
~ 300 K for the 10 /im silicate feature and ~ 150 K for the 20 /im silicate feature for 
example. Thus, 20 fxm silicate features are weak in ULIRGs compared to 10 /xm features 
owing to the dust radial temperature gradient. Ice absorption, on the other hand, should 
only occur in the cooler outer regions of the dust envelope and so may show the same 
relative feature strengths found in the laboratory unless the absorbing covering factor is a 
strong function of the emitting source's wavelength dependent size (buried source radius 
r 0-8(ioiix^)°'^(3~)^ pc where T has been eliminated between the Wein Displacement 
and Stefan-Boltzmann Laws, L is the portion of IR luminosity corrected for extended star 
formation, A runs from 2 /im to the lesser of either the wavelength where the source is no 
longer optically thick or the peak FIR wavelength and scattering is neglected). MIR ice 
absorption is therefore consistent with a portion of the FIR emission being powered by the 
MIR rather than the optical and UV portion of the electromagnetic spectrum. However, it is 
necessary to be cer tain that the reported absorption is indeed owing to ice. As an example. 



Spoon et al.l ((2002) claim ice absorption in 18 galaxies based on 6 //m absorption including 
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IRAS 17208-0014 the third brightest ULIRG. However OH bond bending needn't occur in 
ice alone. The absence of 3 fim absorption in the spe ctra of the dust cori tinuum dominated 
sources Mkn 231 and IRAS 05189-2524 presented by ISpoon et al.l (120021 ). where veiling by 
starlight is unimportant, suggests that an identification of wat er ice may be pr e matu re for 
these sources and refractory materials such as those studied by iGreenberg et al.l (Il995l ) may 
pr ovide a more conv i ncing model. An apparent 6 /xm dip in the spectrum of Mkn 231 plotted 
by IWeedman et al.l (120051 ) similar in amplitude to the 6 .2 //m hump for which they report 
a measurement tends to confirm the ISpoon et al.l (120021 ) claim of 6 /xm absorption in this 
source. The spectrum of Arp 220 ma y be consistent with a contribution of 12.5 /xm libration 
mode absorption, supportive of the Spoon et al. (120021) identification of ice absorptio n in 



this source but, until recently (llmanishi et al 



2006a 



Risahti et al.ll2006l : Ilmanishill2006l ). the 



only well confirmed case of ice absorption among ULIRGs in their sam ple was that of IRAS 
0018 8-085 6 where ice absorption at 3 /xm confirms the identification (llmanishi &: Maloney 
20031 ). See llmanishi et al.l (l2006bl ) for further well confirmed ice identifications. 



To illustrate this we present Fig. 1 which employs a slightly detailed model to explain 
the absorption features obse rved between 2.7 and 12.4 u m. The upper panel reproduces the 
3 fxm spectrum reported by llmanishi fc Maloney! (120031 ) together with an absorption model 
composed of amorphous ice and refractory carbonaceous absorption apphed to a power-law 
continuum, and veiled by some starburst-like emission. The model is then applied in the 
lower panel with some adjustment made for expected reduced refractory absorption at longer 
wavelengths owing to a presumed dust temperature gradient. It has not been possible to fit 
the entire 3 fim absorption without using ice. Thus, the absorption at 6 fim also is partly 
attributed to water ice. With these data, 12.5 /im water ice absorption, while included 
in the model, does not seem to be required so that one may conclude that two features, 
and particularly the degeneracy breaking short wavelength edge of the 3 /im absorption 
contribute to the fairly robust identification of ice in IRAS 00188-0856. We note that the 
need for strong silic ate absorption in this source is confirmed by the clear presence of 18 /im 
silicate absorption (jlmanishi et al.ll2006bl ) 



Thus, ice absorption could play a significant role in ULIRGs, but the evidence is not 
yet as strong as that showing that the MIR emission is often absorbed by refractory dust. 
Evidence for or against FIR ice emission or absorption in ULIRGs is in worse shape still. This 
is primarily owing to a lack of observational data. In Fig. 2 we present archival ISO LWS 
spectra of Mkn 231 (filled diamonds) and IRAS 17208-0014 (filled circles) which show spectral 
structure in the shortest wavelength d etector band SW 1. T his structure is compared to the 
spectrum of the icy source HD 161796 JHoogzaad et al.ll2002h transformed to absorption. For 
Mkn 231, as just noted, the absence of 3 //m ice a bsorption wo u ld lea d to somewhat strained 
geometric arguments to allow acceptance of the ISpoon et al.l (|2002| ) identification of the 6 
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//m absorption, but once made, such arguments would likely allow absorption at even longer 
wavelengths. A complex geometry has alre ady been invoked for this source on a smaller 
physical scale to explain its low X-ray flux (IMalonev fc Revnoldsll2000f). For IRAS 17208- 
0014, the present spectrum would tend to conflrm the Spoon et al. J 20021) ide ntiflcation of 



the 6 /um absorption as do 3 //m spectra ( Risaliti et al. 20061 : Imanishi 2006 ). For either 



source, the presence of crystalline ice implied by the possible 43 fj,m absorption would imply 
a low cosmic ray abundance at least in the region wher e the absorption is prod uced since ion 
bombardment erases crystalline structure in cold ices (IMoore fc Hudsorull992l ). The inset in 
Fig. 2 shows the relationship bet ween the wavelength of the 43 /im peak and ice temperature 
measured by ISmith et al.l (119941 ) which we h ave invoked to shift t he ice spectrum derived 
from HD 161796, an evolved ice forming star (IHoogzaad et al.ll2002l ). beyond the corrections 
made for the source redshifts to better match the spectra. Inversion of these corrections 
would imply temperatures at which amorphous ice may exist and would be preferentially 
formed. Constraints on cosmic ray abundance would need account for the rate at which 
crystalline ice is introduced into the absorbing region either through (substantial) changes in 
luminosity of the source which may allow the annealing of amorphous ice to crystalline ice, 
processing of ice by local heat sources in the region, or the migration of crystalline ice into 
the region. The discovery of sub stantial radial molecular gas motion in the ULIRG IRAS 
08572-1-3915 (iGeballe et al.ll2006l ) suggests the third possibility needs close scrutiny. Since 
the SWl detector of LWS is problematic, these sources must be independently re-observed 
and the potential features conflrmed before fully accepting that ice is present in these sources. 



2.2. Millimeter morphology 



We have examined the LWS SWl spectra of a two other ULIRGS, Arp 220 and Mkn 
273, flnding no feature consistent with crystal line ice absorption. In the cas e of Arp 220, the 
FIR spectrum appears to be optically thick (iGonzalez- Alfonso et al.ll2004l ). There is every 



reason to think ice should be present owing to the source temperature, the clear presence 
of water vapor, as well as support ive MIR observatio ns already mentioned. Here, ice may 
help to resolve an issue raised by ISoifer et al.l (119991 ) with regard to FIR sources that are 
optically thick showing any detectable emission in the MIR owing to overwhelmingly large 
implied MIR opacity. Since ice strongly enhances the FIR opacity compared to refractory 
dust, a given FIR op tical depth c an be achieved with less dust, up to a factor of seven in the 



models computed by lAannestadl (119751 ) . This effect can make the issue raised by ISoifer et al. 



(11999! ) less troubling, tho ugh other approache s such as examining the effects of the dust 



density distribution {e.g. iNenkova et al.ll2002l ) may also need to be explored as well. Ice 



may also have the effect of making millimeter continuum size estimates smaller than what 
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might be expected for refractory dust that is optically thick in the FIR since the opacity of 
ice mantled grains is lower than that of bare grains at millimeter wavelengths. This could 
yield a cuspy morphology in the millimeter where one would usually expect a relatively 
more flattened morphology owing to the cooler outer regions being relatively brighter in the 
millimeter under the assumption of higher refractory dust opacity. 



2.3. FIR ice spectrum 



The infrared labora tory spectrum of crystalline water ice has been well known for some 
time (IBertie et. allll969l ) and differences between the spectrum of crystalline a nd amorphous 



ice h ave been exploited to aid understanding of astronomical sources {e.g. iDartois et al. 
19981 ). However owing to the low thermal conductivity of amorphous ice, preparation of 
laboratory samples, via vapor deposition, that are thick enough to conduct investigations at 
wavelengths longer than 100 ^m has been too difficult until recently. In particular it has been 
unknown whether or not amorphous ice has a fe a ture c orresponding to the broad feature 
centered near 150 /im identified by iBertie et. al Jl969h as a phonon mode in crystalline 
ice. Since nature provides sources rich in amorphous ice, it may be that such a comparison 
could be made through observations of these sources. The spectral coverage of the ISO LWS 
spectrometer extended to ~ 200 /xm, too short to fully cover the bro ad 150 Atm cry s talline 



ice feature. Comparison between the ISO LWS data presented by iDartois et al.l (Il998l ) 



for RAFGL 7009S, a de eply embedded ultraconi pact H II region, and the submillimeter 
observations provided by lMcCutcheon et al.l (119951 ) yield a modified optically thin blackbody 
parameter set of T=27 K, 13=2.1 (where the dust emissivity is proportional to A^^) when 
attempting to the match data near 190 /im. The large value of (3 is suggestive of the presence 
of a feature similar to that of crystalline ice in this amorphous ice dominated source. However, 
full spectroscopic confirr nation is need e d to d ecide if the long wavelength departure of the 
data from the model of IDartois et al.l (119981 ) indicates the existence of an amorphous ice 
feature analogous to the 150 /xm crystalline ice feature. 



Laboratory work by I Curt is et al.l (120051 ) has explored the optical properties of amor- 
phous ice bet ween 15 and 2 00 /xm . There is little question that most of the ice in the thick 
films used by I Curt is et al.l (120051 ) deposited at T < 126 K differs from ice deposited above 
this temperature since the wing of the libration band is clearly shifted so it seems safe to 
conclude that amorphous ice, clearly present in the thinner samples, is also largely present 
in the thicker samples despite issues related to the low t hermal conductiv i ty of amorphous 
ice which hampers the (rapid) growth of thick samples. iGerakines et al.l (120051 ) have also 
had recent success in growing thick samples of amorphous ice using low deposition rates. 
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Schober et al.l (120001 ) take a different approach: annealing samples of high density amor- 



phous ice into the low density amorphous ice found in space . As with the LWS spectroscopy, 
the new amorphous ice optical constants (ICurtis et al.ll2005l ) do not extend to the wavelength 
regime (300 /xm ) where crystalline ice shows absorption whose wavelength dependence is de- 
scribed by (IBertie et. allll969l ). Unlike natural sources, however, where temperature and 
grain size effects make full spectral coverage of a broad feature nearly absolutely necessary 
to ascertain its reality, the controlled lab oratory cond it ions a llow a reasonable extrapolation 
to be considered reliable. The finding of I Curt is et al.l (120051 ) that the imaginary part of the 
refractive indices of amorphous and crystalline ice ar e very similar betwee n 100 and 200 fim 
would appear to be consistent with the suggestion of ISchober et al.l (120001 ) that the order in 
amorphous ice indicated by sharp features in inelastic X-ray scattering data may be owing 
to a hydrogen network that is also present in crystalline ice. If this order is indeed present, 
then a shared spectral structure for the interaction of (FIR) photons with phonons between 
the two types would allow a first order substitution of crystalline ice optical constants for 
the unknown values for amorphous ice beyond 200 /xm. 



2.4. FIR ice emission 



Turning to photometry, and returning to ULIRGs, if we accept the ISpoon et al.l (120021) 



classif icat ion scheme then the spect rum of IRAS 14348-1447 as reported by iGenzel et al. 



(119981 ) or lCharmandaris et al.l (120021 ) must be considered Class II, that is 6 fxm "ice" absorp- 
tion is present and partially filled in by 6.2 /xm polycyclic aromatic hydrocarbon emission. 
This source occupies an extre me position in the dust emissivity (/3) vs. 100 /im optical 
depth (tioo) diagram given by iKlaas et al.l (120011 ) with both high rioo (5) and (3 (2). Since 
the effects of optical depth and emissivity are intertwined, a lower optical depth can imply a 
l arger value of (3. For example, (] = 2.4 and rioo = 0.1 does about as well as the iKlaas et al. 



(I2OOII ) parameters for data longward of 100 fim. However, values of larger th an 2 are 



not generally found except where ice is thought to be presen t (e.g. iLis et al.lll998l) or dust 



temperatures are lower than found in ULIRGs (67 K for the iKlaas et al.l (I2OOII ) fit or 26 K 
at the lower optical dep th and large r tes ted here). In Fig. 3 we compare the fit proposed 
for IRAS 14348-1447 bv lKlaas et al.l (120011) w ith an alternative model. The model shown by 



the dashed line is that of Klaas et al 



(I2OOII). while th at shown by the thick solid line is a 



core-mantle grain model taken from lAannestadl (119751 ) but modified to h ave amorphous ice 
features in the following way: First, opacity data fromlSmith et al.l (119941 ) for crystalline ice 
cooled to 70 K was matched to the lAannestadl (119751 ) curve for Vm/Vc=0.95 by inducing 
a u'^ dependence to approximate the Mie calculation results along wi th the addition of 
an underlying continuum opacity. Second, amorphous ice opacity data (jSmith et al. 1994; 
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deposited at 10 K and warmed to 70 K) was substituted between 30 and 70 /im based on 
this matching after adjustment for the differing sample thicknesses. This opacity curve was 
used to calculate a model with 16 and 71 K co r npone nts. If we take temperature and optical 



depth as free parameters for the iKlaas et al.l (120011 ) fit (2 parameters; normalization fixed 



by IRAS) and two temperatures, two optical depths, and the relative contribution at 90 fim 
of these two components as 5 free parameters for the present model then we may estimate 
how well the two models represent the data if we know something about the errors in the 
data. The statistical errors for the data are thought to be small but the absolute calibra- 
tion error could be large, ~30%. What is important for comparing models is the relative 
accuracy from filter to filter. We have thus assumed 13 % errors for the ISOPHOT dat a to 



give a filter-to- filter accuracy of 20% flKlaas et al.ll2000l ). The fit of iKlaas et al.l (120011 ) for 



data longward of 20 /xm is unacceptable (P(x^,z/)=0.02) under our assumption about the 
errors. A marginally acceptable fit (P{x^,iy)=OA) is found for the amorphous ice model. 
Detailed inspection of the ice fit suggests that the 150 /im ice feature is not accounting for 
all of structure in the data. The thin solid line (Vm/Vc=0.33, T=28.8 K; tqo = 0.2 and 
Vm/Vc=0.95, T=80.5 K, tqq = 0.08 crystalline ice; onion-skin model with warm, T=180 K, 
refractory dust in the interior) shows that this might be possible. However, our method of 
substituting amorphous for crystalline ice is unlikely to work well for the thin mantle used 
here. Given the assumption regarding the data errors made here, it is not possible to make 
a clear claim that ice emission is detected in this source. Assuming larger errors could make 
both low free parameter models acceptable while reducing the errors could make the simpler 
ice model unacceptable while the more complex model cannot be evaluated via reduced x^- 
Clearly, Spitzer MIPS SED mode observations could be definitive since the strong 70 /xm 
curvature in the ice models cannot be produced with standard refractory dust. 

Thus, both astronomical and laboratory far infrared data leave questions about ice 
unanswered. Crystalline ice has a feature near 150 /xm which has yet to be completely 
compared with amorphous ice. Comparison between crystalline and amorphous ice has been 
made near 40 and 60 /im and they are quite distinct whereas they are similar between 100 and 
200 /im. Beyond 200 fim questions remain which are presently answered by extrapolation 
of the laboratory data and hints from the astronomical data. Substituting crystalline ice for 
more common (in space) amorphous ice to consider wavelengths longer than 100 fim could 
add to (systematic) uncertainties in photometric redshift estimation that can not be well 
addressed without further observations or measurements. Regardless of the remaining issues 
concerning the 150 /im feature, for the main purpose of this work it is prudent to consider 
the more extreme consequences of crystalline ice emission on redshift estimation which are 
examined after the next section. 
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3. Ice at high redshift 
3.1. Pair instability SNe enrichment 



Just as interplanetary grains may preserve a record of the sources of ISM enrichment 
occurring at what is presently observed to be z = 0.43, the solid phase of the ISM at high 
redshift may also record the conditions of the earliest enrich ment. Particularly for sources 
that contribute to the reionization of the universe {7 < z < 14; ISpergel et al.ll2006l ). scenarios 
which attribute substantial enrichment to pair instability SNe early on may leave a mark on 
the solid phase of the ISM which could linger through some subsequent processing. A chief 
attribute of pair instability SN enrichment is an early high abundance of oxygen relative to 
other metals that usually deplete to the the solid phase of the ISM. 

For example, a 186 M^r^ progenitor yield s number abundance ratios of 0:Si:C:S:Mg:^^Ni 
= 126:29:15:10:8:1 Jfleger fc WooslevI l2002h which, after formation of CO (29 O removed), 
olivine (16 O and 4 Si removed), silica (50 O removed), and SO2 (20 O removed), leaves 9% 
of the oxygen available to form OH and H2O. For a 100 Mq progenitor, 0:C:Mg=48:9:l with 
little Si, S or ^^Ni leaving 80% of the oxygen available to combine with hydrogen assuming 
the refractory solid state is MgO. In the absence of ^^Ni, the formation of solids should be 
enhanced since it may commence at a higher density in the SNR. These estimates for the 
available oxygen are essentially lower li mits, more complete treatments which con sider the 
possibility that the ejecta are unmixed (INozawa et al.ll2003l : ISchneider et al.ll2004l : allowing 
the formation of graphite for example) would leave an even larger fraction of oxygen available. 
Thus, we might expect grains with large mantle-to-core volume ratios to result from a period 
of pair instability SN enrichment. The continued presence of these ice grains would depend 
on the rate of reprocessing of this initial solid phase in the ISM. 



3.2. Ice emission in a high redshift quasar 



Bertoldi et al. 


2003bh sue- 


of reionization (IFan et al. 



2OO3I ). That the quasar exists indicates that the source has been contributing to reioniza- 



tion for some time and we may take the beginning of reionization as an estimate of it's age 
(~500 Myr) which gives a high average accretion rate for a blackhole mass of 3 x 10^ Mq 
(IWillot et al.ll2003l ). Thus, it seems plausible that pair instability SN enriched material with 
a high ice content is being supplied toward the center of this activity as the blackhole and 
presumably its stellar system grow in mass. On the other hand, ice also forms in the more 
evolved ISM of the Galaxy so an indication of the presence of ice needn't be fully attributed 



- 11 - 



to ice rich regions that might be formed by pair instabihty SNe. 

Fig. 4 shows an indication that the rest frame FIR emission from QSO J1148+5251 
is better explained by ice than refractory dust emission. Three models are compared to 
photometric observation s: 1) a refrac t ory rn odel (dashed line T=55 K; tqq = 0.3; = 2) is 



similar to that given by iBeelen et al.l (120061 ). 2) a simple ice model (dot-dashed line; T=43 
K; Teo = 0.4, Vm/Vc=5.6 amorphous), and 3) a more complex ice model described below 
(solid hne). 

In Fig. 4, it is not possible to formally access the goodness of fit through a reduced 
xl calculation even for the refractory dust model (dashed line) since the number of free 
parameters (normalization, temperature, optical depth, and opacity index) is the same as 
the number of points to be fit. If one were to assume (3 = 2 owing to insufficient time 
since the first enrichment to produce large grains then = 10.7 which strongly rules 
out (P(x, '^)=0.001) this model. That the fit is done by eye probably does not affect this 
assessment but a failure to account for systematic errors could since different groups give 850 
yum photometric estimates using the same data that vary by as much as the refractory dust 
model differs from the present data. By selecting a particular mantle-to-core volume ratio 
(5.6) we may comparably provide a single degree of freedom and obtain = 3.9 (dot-dashed 
line) which would indicate that the model is not fully ruled out by the data (P(x, z^)=0.05) 
and is an improvement on the refractory dust model, again with the proviso that we have 
only considered statistical errors. More elaborate (multi-temperature) refractory dust or ice 
mantled grain models can fit the detections essentially perfectly and are useful to calculate 
as they can indicate what further observations would be decisive. The solid line in Fig. 4 is 
such a model where T=71 K; rgo = 0.4 ice-mantled dust is obscured by T=24 K; tqq = 0.4 
ice-mantled dust. A similar model using j3 = 2 refractory dust (not plotted; T=300 K; 
Teo = 0.1 is obscured by T=1 8 K; rm = 1 dust) e xceeds the plotted 400 /im upper hmit 



arrow in Fig. 4; 0.39 mJy 3(t; iBertoldi et al.ll2003al ) by a factor of 1.9 and is unphysical in a 



number of ways including that the cool component is cooler than the CMB aX z = 6.42. The 
ice models distinguish themselves from the plotted refractory dust model by being fainter 
by at least a factor of 4 at rest wavelength 400 ^m (3 mm observed frame). The present 
upper limit does not distinguish these cases. Such a measurement should be possible using 
CARMA. Determining the composition of the emitting material can also affect estimates of 
the dust mass and thus estimates of the amount of past enrichment (and star formation) in 
the QSO J1148+5251 system. For the present ice model, 7 times as much refractory dust 
would be needed to produce the same FIR opacity so that a lower level of total enrichment 
might be accommodated. 



Finally, at high redshift, where theory suggests that we might most anticipate it owing 
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to expectations about early ISM enrichment, we find evidence for FIR ice emission which, 
while still debatable, seems moderately persuasive and which can be relatively easily checked 
with further observations using existing telescopes. Now we may turn to our main theme. 
We have employed amorphous ice to consider IRAS 14348-1447 and QSO J1148+5251, but 
now we consider crystalline ice both because its effect on redshift estimation is the most 
dramatic and because ice may form at fairly warm temperatures in the SNRs associated 
with pair instability SNe. 



3.3. Ice emission in 450 fim submillimeter sources 



A set of so urces that could be at high redshift and which demonstrate the effects of ice 
are reported by iPope et al.l (120051 ). They are detecte d at 450 fj,m but not at 850 fim. We 
co nsider the 7 s ource s tabulated by iPope et al.l (120051 ) and t wo sources f r om th e original list 
of iBorys et al.l (120031 ) which were recovered above 3.5 a by iPope et al.l (120051 ). These have 
drawn our attenti on because the reported 850 //m upper limits indicate large values of (3 
( IBorys et al.ll2004j : based on original list). Large values of f3 can indicate the presence of ice 
emission. If we generate 3a lower limits to the 450-to-850 iJ,m flux density ratios using the 450 
//m measurements less 3/ a/(2) x CT450 together with 3/ a/(2) x ctsso 850 fim upper limits then 
the lowest lower limit is 8.7 and is only one of two that can be consistent with (3 <2 dust. 
This method of calculating limits can be more conservative than calculating the ratio of the 
450 measurement to the 850 /xm 3(t upper limit (giving 13.4 in this case) but its extension to 
2a should be viewed with caution since detections at the 2.9 a level are customarily reported 
as 3(7 limits. In such situations, when data with comparable signal-to-noise ratios become 
available, the "2(T extension" may be violated much more frequen tly than the name implies 
if many of these near misses lurk in the o riginal data. Howev e r. iPope et al.l (120051 ) report 
no detection from applying the method of ICaillault fc Helfand Jl985h (stacking) to the 850 
/xm data at the positions of the 450 /im sources they detect above 3.5cr. In five cases where 
upper limits are calculated below we substitute the ratio of the 450 /xm 3(j upper limit-to-the 
850 measurement when it gives a more conservative value. 

In Fig. 5 we plot, as a function of redshift, the 450-to-850 ^m flux density ratio calcu- 
lated for a series of SEDs. The dashed line is for a model with T=60 K, 13=2, and T9o=2. On 
this line are plotted at z = 0.9 and z = 0.5 the 3a lower limits (filled circle and diamond) 
for SMM J123603+620942 and J123631+620657. The arrows in the figur e extend to the 
2a lower limits. In these cases the argument given by iBorys et al.l (120031 ) for z < 1 may 
be applied directly, though tenuously, at the 3a level with no assumptions about ice but 
this is not the case for the 7 remaining lower limits. There may be some difficulty however 
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accommodating these two sources at these low redshifts unles s sources with 0.2-2 Jy level 
flux densities are present at 70 fim while iFrayer et al.l (120061 ) find about four ~ 30 mJy 
sources at 70 iJ,m in this region and we are unaware of any IRAS sources in this field. At 
these low redshifts, it is a losing proposition to attempt to increase the dust temperature 
to increase the redshift to avoid the 70 /im constraint in an effort to avoid requiring ice 
emission. When ice is invoked, then these and the remaining lower limits may all be ac- 
commodated but the range of consistent redshifts expands. The remaining curves in Fig. 
5 correspond to increasing ice temperature. The dot-dashed curve uses an SED similar to 
that shown in Fig. 1 (T=30 K), the triple-dot-dashed curve has T=50 K as in Fig. 7., 
and the sohd curve has T=150 K; tqq = 0.04 and Vm/Vc=5.6 for all ice models. For the 
solid curve, ambiguity between amorphous and crystalline ice is physically resolved in favor 
of crystalline ice. The SMM J123603+620942 and J123631+620657 limits are replotted on 
the solid curve giving a redshift limit z < 10.2. This limit is little changed if the maximum 
possible ice temperature (170 K) is used. All the 3a limits can be accommodated at z = 9.5 
but only SMM J123603+620942 and J123631+620657 can also be accommodated there at 
the 2a level. SMM J123657+622033 (filled square) shows this clearly where it is plotted at 
z = 9.6. The 2a level for SMM J133747+621600 (filled clover leaf) also constrains the ice 
temperature to be T^40 K. It is plotted on the T=30 K curve at 2; = 1.9. The remaining 
sources other than SMM J123727+621042 are plotted with redshift upper limits well within 
the epoch of reionization. However, as with sources plotted on the T = 30 K curve, the 
2a limits cannot be accommodated for z ^ 2.5. That 7 of 9 of the 450 //m sources in the 
GOODS-N field require ice to explain the 850 /xm non-detections without recourse to any 
other constraint suggests strongly that the possible effects of ice emission must be considered 
when interpreting submillimeter observations. Ice dramatically changes redshift limits: by a 
factor of 20 for SMM 123631+620657. 

From the other direction, for the so urces with 850 nm detections and 450 /im upper 
limits that also have estimated redshifts (IBorys et al.ll2004j ) we may also begin to constrain 
the presence of ice. 14 upper limits are plotted in Fig. 5 as open circles for sources where 
a redshift estimate is given. The open diamonds are also upper limits plotted (but withou t 
arrows to avoid crowding) at the lower limits to their redshift given by iBorys et al.l (120041 ). 
For redshift estimates with z < 3 the upper limits are constraining on the likelihood that 
optically thin ice contributes to the FIR opacity: at the 2a level (arrow heads) little T>30 
K optically thin ice seems to be needed. At the 3a level 4 of 11 sources could be consistent 
with a contribution of optically thin ice emission to the source SEDs. One source with a 
redshift lower limit {z > 2.7; SMM J123652+621225) needs to be optical ly thick and cool 
(T ~40 K) or at a higher redshift as already noted by iBorys et al.l (120041 ). The remaining 
sources with with redshift estimates above 3 or redshift lower limits are not constraining 
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on the presence of ice. The present considerations suggest that the 450 and {z < 3) 850 
fim detected sources differ physically in that the former have optically thin ice emission that 
dominates their FIR opacity while the latter may be optically thick at 850/(1+2;) //m or that 
refractory dust dominates their FIR opacity or both. The average 450-to-850 /xm ratio for 
the 850 /im source c alculated fro r n the stacked 450 /xm and the average of the 850 /im flux 
densities reported by lPope et al.l (120051 ) is 0.5 which we have plotted as a small open square 
a.t z = 4.25. The plotted redshift would be consistent with T=30 optically thick material 
or T=20 optically thin ice or (3 = 2 dust. If it is true that 850 /im sources are on average 
fainter at 450 /im than at 850 fim then it seems at though many of the 850 /im sources 
must lie a.t z > 4 and ice may be needed to make them detectable at 850 /xm with ULIRG 
luminosities (see end of Sec. 4.2). We do worry, however, that the stacked 450 /im flux is 
underestimated owing to the smaller 450 /tm beam or that the data were normalized prior 
to stacking making the reported units ambiguous. QSO Jl 148-5251 is also plotted (small 
filled circle). As discussed above, it seems to require ice emission. 

Just as detected 850 /im sources must be at least ultraluminous, so also must the 
detected 450 /im sources be hyperluminous (L> 10^^ L©). In Fig. 6 we plot the models 
shown in Fig. 5 giving their 450 /im brightness as a function of redshift under the following 
assumptions about luminosity: The models for refractory dust, and ice at T=30 and 50 K are 
ten times more luminous than ULIRG IRAS 00188-0856 (Fig. 1), the model for T=150 K ice 
is 10 times more luminous still, and we have also plotted the complex amorphous ice model 
for QSO J1148-I-5251 at it's observed luminosity (short dashed line and filled circle; ^ 10^^ 



Lq] iBeelen et al.ll2006l ). The brightnesses of the 450 /im sources given by lPope et al.l (120051 ) 



lie in the range 77-to-291 mJy so that the plotted models would need to be scaled up by a 
factor of a few to cover this range fully unless the sources are at 2; < 0.2. For example, one of 
the brighter sources (SMM J123631-I-620657; 263 mJy) would have a luminosity of 6x10^^ 
Lq for the T = 150 K ice model if located at z = 9.5, a factor of three higher than the most 
luminous hyperluminous source listed by IVerma et al.l (l2002l ). It's luminosity would be 10 
times lower for the T = 30 K ice model located at z = 2, near the middle of the distribution 
of luminosities in their list. From Fig. 5 we estimate that the 450 /im sources would populate 
the redshift ranges 1-3 and 5-10, with the latter range having a larger co moving volume by 
a factor of 1.7. If ice is only important close to reionization so that the sources were all at 
9.2< z < 10.2 (see solid line in Fig. 6) then the luminosities require about a tenth of the 
stellar mass of M87 per source to be involved in star formation. The sp ace density would be 
2.5x10"^ Mpc~^, comparable to the local density of clusters of galaxies (iBahcall et al.ll2003l ) 
and the luminosity density woul d be ~ 5 x 10^ Lf;^/Mp c^ (comoving), about a factor of 100 
higher than the local FIR value (ISaunders et al.lll990l ). This would then imply a minimum 
of 10^ recombinations per proton if the 450 /im luminosity traces a factor of ten smaller 
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ionizing luminosity: about 10^ more recombinations than for the hfetime of a typical H II 
region which seems too high for reionization unless dust competes efficiently for reionizing 
photons. 

On the other hand, there are six 850 /im sources with flux densities above 20 mJy listed 
by Pope et al. (2005) which could be close to hyperluminous (see Fig. 7) so that we need 
only postulate a similar population which is (thus far) invisible a 850 //m owing to the effects 
of ice and thus expect to flnd these sources at a range of redshifts. 

There is very clear evidence of ice rather than dust providing the FIR opacity in 7 of the 
9 450 /im sources considered here. But, the redshifts may or may not be high and additional 
photometry is needed to constrain the redshifts. 



HDF 850.1 and ice 



4.1. 



The search for counterparts at other wavelengths to the br ightest 850 fim source in the 
Hubble Deep Field has been arduous but iDunlop et al. J2004h make a convincing case for 
detections in the radio and near infrared (NIR) and summarize submillimeter observations. 
Modeling fl Aretxaga et al.ll2003l ) of these observations leads them to conclude that the red- 
shift of HDF 850.1 lies in the range 4.1±0.5. Since the SED of a high q ULIRG (IRAS 
08572+3915) is included in the 20 templates in the model input, the effects of the change 
in the radio-FIR relation for ULIRGs compared to lower luminosity galaxies noted in the 
introduction should be partly represented in their estimate. Similarly, the suggestion of ice 
emission in the SED of IRAS 14348-1447 (Fig. 3) could also influence the estimate if the 
photometric data were used directly rather than fitted with a modified blackbody. In gen- 
eral, reliance on modified blackbodies will give lower estimated redshifts for cooler sources 
when more than one point on the Ray leigh- Jeans curve is available, as is the case for HDF 
850.1, but the addition of an ice emission feature can increase the estimated redshift even 
for a cool source. 



As one example, the bow-tie shown in Fig. 1 on the FIR continuum corresponds to the 
range of reds hift one would calculate for HDF 850. 1 based on two estimates of its 850 fim 
flu x densitv (jDunlop et al.l |2004J : IWang et al.l |2004| ) and the 1.3 mm flux density reported 
' 1999h. That is 5.4 < z <5.6. It should be noted, even before describing 



by iDownes et al.l ( 



the model upon which this is based, that a detailed comparison between the SED shown in 
Fig. 1 and the optical and NIR limits and detections for H DF 850.1 is no t poss ible owing 
to insufficient wavelength coverage for IRAS 00188-0856. See iDunlop et al.l (120041 ) fig. 6 for 
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a comparison with data from two ULIRGs and a luminous infrared galaxy calculated for 
a range of redshifts. However, the effects of far infrared ice emission can be demonstrated 
with this example. The model, the sh ort-dashed line in the lower panel of Fig. 1, is based 
on curves given by Aannestad (jl975 ) which are the results of calculations of the optical 
properties of ice mantled crystalline silicate grains. And, it is crystalline water ice in the 
mantles. We chose the mantle-to-core volume ratio = 5.6 curve to most closely match the 
ratio of 3 //m ice to silicate optical depth used in the model shown by the solid lines in 
Fig. 1. For a number of reasons this choice is probably not fully constrained by this ratio: 
the optical depth to which the silicate absorption is sensitive may be larger than for ice; 
an odd shift in the wavelength of the silicate feature may indicate a pyroxene composition 
with a different feature- to-continuum optical depth ratio than usual; the crystalline silicate 
in the I Aannestad! (119751 ) calculation is somewhat dissimilar to astrophysical silicate models. 
However the first two effects may be compensating so that as a rough abundance analysis, 
the expectation that ice would dominate the FIR opacity in this source seems well founded. 
This condition, and the need for the ice emission to be moderately optically thin are really 
the only two requirements for ice to strongly affect redshift estimation. The optical depth of 
the model of far infrared emission shown in the lower panel of Fig. 1 is set to correspond to 
the 3 fim ice optical depth used in models shown by the solid line in the upper panel (where 
ice and refractory material contribute equally to the absorption optical depth at 3.1 /im). 
The temperature of the far infrared emission model is set to roughly reproduce the IRAS 60 
and 100 fim measurements. This temperature also accentuates the effects of the broad 150 
/im feature over which the HDF 850.1 data are displayed. The range of redshifts needed to 
place HDF 850.1 on this model falls outside of the range estimated by iDunlop et al.l (120041 ) 
and gives a first indication that the effects of ice on photometric redshift estimation (as 
opposed to just limits) can be significant. 



4.2. z = 13 



As noted in the introduction, using local ULIRGs to understand submillimeter sources is 
a sensible approach if their properties are not strongly affected by cosmic evolution. However, 
at higher redshift, this proviso seems fairly likely to break down. Many local ULIRGs 
appear to be produced in the mergers of gas rich spiral galaxies and it has been argued 
( iMihos fc HernquistI 119961 ) that the dynamical stability provided by bulges may preserve 
gas reservoirs for the final merger stages associated with some ULIRGs. Owing to the 
finite lifetime of the universe, there are earlier periods when the precursor systems may not 
resemble present day precursors i n their dynamica l state or star-to-gas mass ratio. It is 
not at all clear that seed SMBHs (lOsterbrockl Il993l ) would be present at sufficient mass to 
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allow accretion to b e an energetically imp ortant or dominant power source as suspected for 
local ULIRGs {e.g. Ilmanishi et al.ll2006al ). Thus, in what follows we assume star formation 
as the main power source at higher redshift. A sketch of how crystalline ice might exist 
in a primordial galaxy environment is also given to explore the plausibility of a very large 
effect of ice on redshift estimation. It is to be noted that giving up a central power source 
such as would be provided by an AGN would usually mean that ice mantled grains would 
contribute little to the FIR emission since most of the energy that powers FIR emission would 
be absorbed in reactive photodissociation regions (PDRs) rather than in molecular clouds 
where ice is thought to exist unless ar i improbable s cenari o involving an overabundance 
of protostars such as that propo sed by iRoussel et al.l (120031 ) but previously discounted by 
Dudley &: Wynn- Williams! (119971 ) is invoked. Below we argue for such an overabundance. 



In Fig. 7 we sh ow the observations of HDF 850.1 at z = 12.6 This is 18 confidence 



intervals beyond the iDunlop et al.l (120041 ) estimate. So far as we can tell, neither the new 
Z bar id upper limit reporte d here, nor the released Spitzer IRAC and MIPS GOODS im- 
ages (IDickinson et al.l |2006| ) rule out either the model shown in Fig. 7 or the analysis of 
Dunlop et al.l (l2004l ) so long as a reddish ULIRG (see iDudleyl (Il999) for a discus s ion o f 
12-to-60 / xm flux density ratio s in ULIRGs) is used as a basis for the iDunlop et al.l (120041 ) 
analysis. IDunlop et al.l (120041 ) argue that the lensing source (solid diamonds in Fig. 7) is 
evolved, in which case much of the tentative observed frame 24 flux density could be 
owing to HDF 850.1 but little of the 8 /xm flux density. At ^ = 4.1, Mkn 231 and 273 could 
be too blue to account for this level of emission whereas Arp 220 or IRAS 17208-0014 would 
not violate this datum if treated as an upper limit. On the other hand, at 2; = 1.1 the chances 
of finding a younger elliptical galaxy to act as a lens are increased so that some or all of the 
observe d 24 nm flux density could be owing to the lensing galaxy. This would still not rule 
out the IDunlop et al.l (120041 ) redshift estimate for HDF 850.1 but would accommodate the 
present estimate under the assumption of the plotted starburst model. The shaded region in 
Fig. 7 shows a possib le range for the spe ctral shape of the lens source emission. The source 
labeled lb in fig. 8 of IWang et al.l (120041 ) may also contribute to the 24 /xm emission though 
only slightly to the IRAC d ata. If it is similar to NGC 253, and at a redshift of = 1.76 
( iFernandez-Soto et al.lll999l ) then up to ~30% of the 24 yum fiux density might be attributed 
to this source but again negligibly at 850 //m. 

The model for HDF 850.1 shown in Fig. 7 (thick gray line) is similar to that shown 
in Fig. 1 in that it has the same mantle-to-core volume ratio and optical depth however 
the temperature is higher as it must be owing to the higher CMB temperature at 2; = 12.6. 
Whereas in Fig. 1, the 150 /im ice feature spectral structure infiuences the redshift placement 
of the HDF 850.1 data, in Fig. 7, shorter wavelength spectral structure also owing to ice is 
infiuential. Interestingly, assuming that there is a foreground neutral intergalactic medium 
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as a result of incomplete reionization (jSpergel et al.ll2006l ). one would expect Gunn- Peterson 
saturation or a Lya blackout to occur in the H band for z = 12.6 the shortest wavelength 
band for which HDF 850.1 is detected. This is represented by the change from black to gray in 
the thin line using a Heaviside function as multiplier. Here th e thin line is a young starburst 
spectrum generated using Starburst99 Jleitherer et al.lll999h. It is notab le that taking the 
onset in redshift of the Gunn- Peterson trough to be ~6 (IFan et al.ll2004l ) absorption at the 
Lyman limit at 2; < 12.6 easily overlaps Lya absorption between z = 6 and z = 9.2. One 
expects pretty thorough opacity shortward of Lya. 

While the choices of parameters of the starburst model are meant only to be suggestive, 
for the model shown in Fig. 7, the free-free emission expected based on the rate of produc- 
tion of ionizing pho tons would explain about 20% of the observed radio emission (following 



Scoville et al 



I991I ) while retaining the starburst model as it is and correcting for sufficient 



extinction to account for the model infrared emission would over produce the radio emission 
so that a lower contribution of stars which produce ionizing photons to the observed frame 
band flux might provide a more likely scenario. Alternatively, invoking a chemical trig- 
ger for star formation discussed below, it is conceivable that an object could be caught in 
a phase when a sufficient number of stars are in an ultracompact H II region configuration 
that, given a sufficiently high ionization parameter, higher than usually required if the gas- 
to- dus^jii^ssj;atioisJiigh owing to low enrichment, free- free emission would be suppressed 



.e.g. 



Bottorff et al.lll998h . 



As noted in the introduction, sources at high redshift should eventually drop out of 850 
/im surveys. In Fig. 8 we plot, for the same far infrared luminosity, the expected 850 /im 
brightness as a function of redshift for an Arp 220-like SED and the two distributions shown 
in Figs. 1 and 7. The 30 K model is truncated at 2; = 11 since the cosmic background 
radiation temperature exceeds the model temperature above this redshift. As can be seen, 
the models with ice select for high redshift sources assuming such sources exist. For the 30 
K model the brightest (non-local) source would be found near a redshift of 5 and for the 50 
K model near a redshift of 13. Should ice play the role proposed here in a moderate fraction 
of the detected 850 /im sources, it would be of little surprise that the brightest source in the 
HDF (850.1) would be among them. 



4.3. Why star formation? 

Our choice of presenting a model that uses star formation as a power source may seem 
strange given that PDRs would be inefficient at producing ice emission which we require to 
allow, for example, HDF 850.1 to be detected at z = 13. One phase of star formation that 
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does allow MIR radiation to heat ice is the protostar phase. This phase is brief compared to 
O star life times and thus would not typically be expected to provide the dominant source 
of FIR radiation in a galaxy. Protogalaxies, however, are by definition required to be young, 
so that while the argument concerning the relative brevity of the protostar phase cannot be 
ignored, it may be overcome if the star formation rate is increasing rapidly with time. And, 
it is just such sources that one might expect to be observed as they populate the top of the 
infrared luminosity function since they would be the most efficient at converting radiation 
to FIR wavelengths. 

A rapidly increasing star formation rate might be expected in a situation where the gas 
cooling efficiency is increasing owing to ongoing initial enrichment. A positive feedback can 
lead to an exponentially increasing system response. The number at the right hand of Fig. 
8 shows the 850 brightness (in mJy) of the three SEP sh apes considered here if we adopt 
the lensing amplification proposed by iDunlop et al.l (120041 ) for HDF 850.1 for a source with 
the FIR lumi n osity of IRAS 00188-0856. The increase in luminosity required to match the 
Hughes et al.l (Il998l ) estimate for HDF 850.1 would be a factor of 2.3 for the dashed curve 
taken at a redshift of 13. This would then correspond to ~ 4 x 10^^ L0 {e.g. Fig. 7) or 
10^ O stars and 1.5x10^^ GO stars for a Salpeter initial mass function (m~^'^^), or roughly a 
galaxy's worth of star formation. This then would be the limit for the continued exponential 
increase in the star formation rate: the amount of available enrichable gas. Night et al. 
(2005) find 4 objects with 10^^ M0 in stars at 2; = 6 in their largest cosmological simulation 
(142 Mpc; large enough to contain ~1 ULIRG at 2; = 0). Thus, there is some suggestion that 
relatively massive objects (or perhaps regions that will become objects) are forming their 
stars a few dynamical timescales earlier than this. Notwithstanding the existence of QSO 
J1148+5251 and other high redshift quasars, it seems to us that there is a nee d to account for 
massive sources at intermediate redshift which last formed stars near z = 13 (IStockton et al. 



2OO4J : iJimenea I2OOOI ) as well as the apparent mono tonic increase in the comoving rate of 
gamma-ray bursts out to z = 7 (ILe fc Dermerl 120071 ). 



However, in the low-to-moderate redshift universe, we would usually take ice emission as 
an indication of a buried compact power source at least at ULIRG luminosities, and, should 
blackhole growth commonly reach ~ 10^ Mq by z = 13, as must happen at least occasionally 
given the existence of QSO J1148-f5251, then our low-to-moderate redshift explanation could 
also apply during reionization. 
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5. Summary and prospects 

Far infrared ice emission may be crucially important to understanding the far infrared 
emission of sources at high redshift and may strongly affect the range of redshifts estimated 
for submillimeter sources, even placing them at reionization redshifts. Observations of in- 
frared ice emission are definitive in Galactic protostars and the winds of evolved stars but 
are scant and not fully persuasive for the local analogs of the submillimeter sources, the 
ULIRGs. This, despite mounting evidence that ice is involved in absorbing much of the 
energy that is finally emitted in the far infrared in a substantial fraction of ULIRGs. 

We find stronger indications that ice emission is important at high redshift . The data 



for QSO J1148+5251 as well as the 450 fim sources cataloged by iPope et al.l (120051 ) are 
both strongly suggestive. We feel that further investigation along these lines could be quite 
fruitful. 

We've considered two models where the presence of far-infrared emission features owing 
to ice in the S EP of HDF 85 .1 wo uld place it at redshifts larger than the range estimated 



without ice by iDunlop et al.l (120041 ). In the first, the use of crystalline ice at 150 fim as a 
substitute for amorphous ice is probably justified, subject to further laboratory investigations 
of amorphous ice. In the second, consideration of a first chemical enrichment mechanism 
suggests that crystalline ice could be present at early times while energetic and timescale 
considerations may be compatible with protostars in protogalaxies producing the sort of 
features we consider here. In the second model, a Lyman a blackout would be consistent 
with the NIR photometry for HDF 850.1. Howev er, we do not prefer this model and its 



higher redshift over the redshift range estimated by lDunlop et al.l (120041 ) or z ~ 5.5 from our 
first model. What we find to be important is that far infrared ice emission could have such a 
dramatic effect on the estimated redshift. It should be noted that redshifts estimated from 
the radio-FIR relation could also turn out to be lower when ice emission is invoked since 
a large submillimeter-to-radio flux density ratio is easier to produce with high f3 at some 
redshifts. So, the main conclusion of this work is that ice broadens the range of redshifts to 
be considered for a given submillimeter source and that ice is quite likely to be important 
at high redshift. 

Some astronomical observations which could aid progress on open issues are: 



• Continuum observations at 3 mm of QSO Jl 148+5251 may confirm ice emission in this 
source. 

• The Atacama Large Millimeter Array may be able to isolate the transition from dust to 
free- free/synchrotron emission for sources of this type and allow more reliable redshift 
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estimates. 

• Investigating the role of ice emission locally could be taken up using SOFIA or Spitzer. 

• Very low resolution H-band spectroscopy might test the possibility of a Lyman a 
blackout for HDF 850.1. 

• Molecular studies which target H2O and OH maser emission may give definite red- 
shifts using a new generation of radio telescope arrays even if the continuum is below 
detection limits. 

Colleagues have commented helpfully at various stages of this work. Some are: S. 
Chubb, E. Dwek, J. Fischer, P. Gcrakincs, A. Li, M. Moore, L. J. Rickard, E. Shettle, 
and H. Smith and an anonymous referee. Research in infrared astronomy and particularly 
the astrophysics of celestial backgrounds at the Naval Research Laboratory is supported by 
the Office of Naval Research (USA). This pubhcation makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and the National Science Foundation. 
Use is also made of the USNOFS Image and Catalog Archive operated by the United States 
Naval Observatory, the archive of results from the ESA sponsored Infrared Space Observatory 
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Fig. 1. — Upper Panel: Data from llmanishi fc Maloneyl (120031 ) is reproduced for IRAS 
00188-0856 together with a model (solid line) subsequently employed in the lower panel. 
The model consists of an a = 2 (flux density f^y oc iy^°') power-law continuum to which 
absorption is applied together with an unabsorbed component which is a combination of an 
archival spectrum of M 82 (ISO SWAA11600319) and a blackbody of temperature 2500 K to 
represent additional stellar photospheric emission (displayed in the lower panel long-dashed 
line). The two absorption componen ts used are: (1) a volatile admixture (H20:CO=2:l) 
of amorphous ice mea sured at 10 K (IHagen et al.lll983l ) and (2) the refractory EURECA- 



B curve presented by iGreenberg et al. 



( I995I ). The effects of absorption by the refractory 
material are shown with the dot-dashed line and those of the combined volatile and refractory 
absorption with the dashed line. The volatile and refractory components have approximately 
equal contributions to the absorption at 3 /im but the refractory component contributes little 
at the short wavelength edge suggesting that a volatile component is required. Lower Panel: 
The same model (solid line), now adjusted to lower optical depth in the refractory component 
to better match the 6 /im absorption and with a smaller additive M 82 component is applied 
to longer wavele ngth data. Here, the effect s of sil icate absorption have also been included 
based on fig. 5 of iDudley &: Wynn- Williams! (119971 ) but oddly shifted to the red by about 0.6 
lim at 10 /xm, perhaps owing to differing composition. The inset shows the effects of adding 
the spectrum of M 82 as in the upper panel and without this component in the model at 
all (bounding dotted lines). The spectral resolution of the model matches that of the data 
in the inset. The short-dashed line is an emission model described more fully in the text 
(Sec. 4.1). Data sources fo r IRAS 00188 - 0856 i n descending order i n the legend a r e from 



Imanishi fc MalonevI (120031): Monet et all (jl998|); ICutri et al.l (l2003f ): iMoshir et al.l (Il990f ): 



and Charmandaris et al 



(I20021). R and K ' data are in rough agreement with the photometry 



given by iKim. Veilleux fc Sanders! (120021 ). 



Fig. 2. — Highly processed data products from the ISO Archive (hpdp_05100540_2.tar and 
hpdp_65000608_2.tar with minimal ISAP post processing) are compared with a 43 fim ice 
absorption feature derived from the spectrum of HD 161796 (also from the ISO archive: 
SWAA52100537). The diamonds and circles are data for Mkn 231 and IRAS 17208-0014 
respectively. The solid lines are derived from the ISO SWS spectrum of HD 161796 using 
a sing le temperature (140 K) to estimate emissivity after the manner of !Roche fc Aitken 
(!l984j ) and then converted to absorption. The inset is described in the text. 
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Fig. 3. — Observational data for IRAS 14348-14 47 (filled circles) are plotted at their rest 
wavelength in Jy as reported by lKlaas et al.l (120011 ) . Horizontal bars indicate filter passbands. 
The solid line shows an optically thin (in the FIR) ice emission model whose (3 parameter 
(100 and 300-900/im) is ~ 2.3 as compared to ~ 1.9 for a bare (refractory) core. Adjustment 
of the core to smaller /3 would result in a smaller f3 for the mantled grain as well. In fitting 
the model data at 95 and 450 /xm were matched and various mantle-to-core volume ratios 
were compared with photornetry b etween these wavelengths. The dashed line reproduces the 
fit proposed by iKlaas et al.l (120011 ) which is quite optically thick in the FIR and has /?=2. 



Fig. 4. — Observational data for QSO J115 4-I-5251 (filled c i rcles a nd arrow) are plotted at 



2003ah (3 mm): iBertoldi et al. 
Beelen et al.l (|2006|) (350 /xm). 



their re st waveleng th in mJy as reporte d by iBertoldi et al.l ( 
J2003bh (1.2 mm); iRobson et all J2004h (850 and 450 ^m) 
The lines represent emission models. The dashed and dot-dashed are single temperature 
models using refractory and ice mantled dust respectively while the solid line two temperature 
ice model with colder ice obscuring warmer ice. 



Fig. 5. — The observed frame 450-to-850 flux density ratio calculated for a series of emission 
models is plotted as a function of redshift together with 9 observed lower limits (large filled 
symbols), 34 observed upper hmits (medium open symbols), and two detections (smallest 
symbols). Individual 450 fim source and detection identifications are given in the legend. 
The models, also in the legend, are more fully described in the text. 



Fig. 6. — Brightness in mJy at the observed wavelength 450 /xm is plotted as a function of 
redshift for a set of models. Four of the models are carried over from Fig. 5: The solid line 
represents 150 K ice, the triple dot dashed line, 50 K ice, the dot dashed line 30 K ice and 
the long dashed line 60 K refractory dust. A fifth model (short dashed line) is taken from 
the solid line in Fig. 4: 71 K ice is obscured by 24 K ice. Infrared luminosities for the models 
are constant as a function of redshift and are 2 x 10^^ L©: long-dashed, dot-dashed, and 
triple-dot-dashed lines, 2 x 10^^ Lq: solid line, and 1 x 10^^ L©: short-dashed line. We remark 
that just as standard dust shows roughly redshift independent brightness at 850 fim, so too 
does T=50 K ice (triple dot dashed line) at 450 /im. The filled circle is QSO J1148-I-5251. 
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Fig. 7. — Observational data for HDF 850.1 (filled circles) are plotted for a rest wavelength 
and luminosity scale corresponding to a redsliift of 12.6 and accounting for the effects of 



gravitational lensing by reducing the apparent luminosi ty by a factor o 
for the Z band limit, are reported or su mmarized in 



Dun 



limit is measured from GOODS ACS data (iGiavalisco et al 



3. T he data, except 
OP et all fl200l . The Z band 



20041 ) by subtracting the image 



of the lensing galaxy from itself after a rotation of 180°. The lir nit is 2 a drawri from the 
weight image for an O'.'S aperture centered at the IRAM position (IDownes et al.lll999l ) and 
corresponds to Z — K' > 3.5 (Vega). Data for the lensing source (filled diamonds) should 
be read relative to the HDF 850.1 data (filled circles) as the scalings in wa v elengt h and 
luminosity do not apply. Here the I — K' photometry are from iDunlop et al.l (12004r) while 



the 3 .6-24 /xm points are measured from GOODS IRAC and MIPS images (IDickinson et al. 



20061 ). The IRAC data are 6'.'0 diameter aperture photometry subtracting an equal area 



outer annulus. The 24 /im datum is obtained through radial profile fitting after subtraction 
of a confusing source and is p erhaps a 2.5cr measurement of 20 /xJy. The dashed hne is scaled 
{z = 1.1 wrt observed frame: iFernandez-Soto et al.lll999r) data for the elliptical galaxy NGC 



1399 (IPersson et al. 



1979 : Athev et al. 2002[) and the t riple-dot-dashed line is for the galaxy 



NGC 3928 (lHuchralll977l : 2MASS; IXilouris et al.ll2004l ) with the shaded region between them 
showing the age de pendent range of mid-infrared flux suggested by models for elliptical 
galaxies reported by lTemi et al.l (120051 ). The models shown as solid lines are described more 
fully in the text. 



Fig. 8. — The 850 yum brightness as a function of redshift is shown for three calculated 
SEDs. The solid line corresponds to the FIR SED shown in Fig. 1. The long-dashed 
line corresponds to an optically thin T = 45, (3 = 1.2 modified blackbody similar to the 



Scoville et al.l (jl99lh fit for Arp 220. The short-dashed line corresponds to the FIR SED 



shown in Fig. 7 all normalized to the FIR luminosity of IRAS 00188-0856. The arrows for 
th e sohd, dot-da s hed, a nd dashed lines indicate redshifts for HDF 850.1 estimated using Fig. 
1, iDunlop et al.l (120041 ). and Fig. 7 respectively. The numbers at the right are a conversion 
of the vertical scale for the gravitational lens amplification factor of three estimated for HDF 
850.1. 
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